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The active-sterile neutrino conversion is studied for neutrino propagating in the axial 
potential generated by magnetised electron plasma in the supernova medium. We consider 
the effect of random magnetic field B rm3 on the average neutrino conversion probability. We 
obtained the constraint on Am 2 and sin 2 20 for different strength of the random magnetic 



oo 



3 



X 



fields, by considering the positive definiteness of the average neutrino conversion probability 
inside the supernova core. Our calculation shows that, B rms <^ 0.63 x 10 14 Gauss is preferable 
so that small values of Am 2 can not be excluded. 
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CN . I. INTRODUCTION 

The COBE measurement of the cosmic microwave background temperature anisotropies on large scale 
| Q, hints for the existence of a hot dark matter (HDM) component of 30% of the total mass density 0. 

Also the recent results of solar neutrino, atmospheric neutrino and the LSND result, hints for the neutrino 
oscillation. All these can not be explained with the three generation neutrino flavor mixing scheme. Thus 
the speculation is that, there is a fourth generation low mass sterile neutrino v s which will mix with the 
standard model neutrinos. This neutrino has to be sterile with respect to electroweak interactions so that, 
it can not be detected directly through experimental searches. This postulate of a fourth generation light- 
O 1 sterile neutrino can explain simultaneously the dark matter, solar neutrino and atmospheric problems 



The most stringent constraints for the neutrino mass matrix including a sterile neutrino species are obtained 
from the nucleosynthesis bound on the maximum number of extra neutrino species that can reach thermal 
equilibrium before nucleosynthesis and change the primordially produced helium abundance [Q . 

Neutrino propagation in the magnetised medium has interesting consequences in the astrophysical and 
cosmological scenarios. Large scale magnetic fields in the early universe hot plasma and in the core of 
the supernova can effect the neutrino conversion It has been shown that random magnetic fields can 
strongly influence neutrino conversion rates and this could have important implications, especially in the 
case of conversion involving a light sterile neutrino Pljioll . The effect of active-sterile neutrino conversions in 
a supernova has also been discussed, both in the case where no magnetic field is present, as well as in the 
presence of random magnetic field as large as 10 16 Gauss. 

The dispersion relation for neutrino propagating in a magnetised medium is different from the vacuum one. 
In the magnetised medium neutrino acquires an axial potential which is proportional to the scalar product 
of the neutrino momentum and the magnetic field vector (k.B). The effect of axial potential on neutrino 
propagation in media with regular and/or random magnetic fields are considered in the literature JTo| |l3|| . 
For random magnetic fields case Semikoz and Valle [l3]| have show that, the neutrino conversion is aperiodic 
in nature. The effect of axial potential on active-sterile neutrino conversion has also been considered in 
supernova and sun (l(],|l4[]. It was shown by Kusenko and Segre recently that, the polarisation effects of the 
supernova medium lead to the explanation of the birth velocity of pulsars |15fl . 
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In this paper we have considered the effect of axial potential on neutrino propagation in the supernova 
medium in the presence of random magnetic fields. We calculate the average neutrino conversion probability 
V, for the process v e — > v s using two different approximations (method-I and II). method-I shows that, for 
both r ii > ujf and Tii < ujf (T\\ is the longitudinal damping parameter) one can obtain finite conversion 
probability for v e — > v s . On the other hand for the positive definiteness of the conversion probability in 
method-II, only L|| < u>f is possible. Using this positive definiteness condition of the average conversion 
probability inside a supernova core, we have showed that B rms < 0.63 x 10 14 Gauss is preferred, because 
large values of B rms exclude small values of Am 2 , which are not supported by present experiments. 

The paper is organised as follows: in Sec. II we calculate the average neutrino conversion probability using 
two different methods in the presence of a randomly fluctuating magnetic field. The comparison of both 
the methods are undertaken in Sec. Ill and necessary condition for the positive definiteness of the average 
conversion probability for the process v e — > v s is obtained in the second method. Using the necessary 
condition for active-sterile neutrino conversion inside a supernova core, we put constraint on the parameters 
Am 2 and sin 2 29 for different values of the random magnetic field in Sec. IV and a short conclusion is drawn 
in Sec. V. 



II. NEUTRINO PROPAGATION IN THE MAGNETISED MEDIUM 



A. method - I 



The evolution equation for a system of two neutrinos v a and v b , where v a is the active one and v b is 
active/sterile one is given by 



i- 



d ( v a \ ( H aa (t) H ab (t) \fv a 



dt\v b J \H ba (t) H bb (t) J\ v b 



(1) 



where the quantity H is in general the potential for the neutrino in the medium which we will discuss soon. 
Let us define the functions R = Re({v*v b )) and / = Im((is*is b }). Then using these in Eq.(Q) we obtain 

R(t) = -H d (t)I(t), I(t) = H ab (t) (2P(i) - 1) + H d R(t), (2) 

and 

P(t) = -2H ab {t)I(t) (3) 

where the function P(t) is the neutrino conversion probability P Ua ^ Ub (t) and Hd = H aa {t) — H bb {t) and dot 
on the top corresponds to derivative with respect to t. Using the Eq.(||) in Eq.(||) we can write 

P(t) = -2H ab (t) f i(h)dh 
Jo 

= -2H ab (t) [ [H ah {t x ){2P{t x ) - 1) + H d {ti)R{tx)\ dh. (4) 
Jo 

Let us consider the neutrino propagation in the medium in the presence of a magnetic field. Then we have 

H d (t) =V- Acos26 + Va X i a i and H ab (t) = (iB ± (t). (5) 
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The quantity V is the difference of neutrino vector potential for v a — ► v b , A = (m\ — m\)/2E = Am 2 /2E and 
E is the neutrino energy and is the neutrino mixing angle. The axial vector potential V ax iai = Me//k.B(i)/fc 
is generated by the mean axial vector current of charged leptons in an external magnetic field |Tl|- |l3| , |l6| . 
For fluctuation in the magnetic field we can write B(t) = Bq + B'(t), where Bq is the constant background 
field and B'(t) is the random fluctuation over it. Then we have 

H d (t) = H d (0) + H' d (t) 

= (V-A + fi eff B {]0 )+fi eff B{ l (t), (6) 

and 

H ab (t) = H ab (0) + H' ab (t) 

= fiB ± o+i^B' ± (t). (7) 

For the neutrino conversion length greater than the domain size i.e. l con v » Lq (where i con „ ~ l/Ww 
and T\y is the weak interaction rate), a neutrino will cross many magnetic field domains before it flips its 
helicity. Thus the neutrino will experience an average field before it flips its hclicity. So one can average the 
propagation equation (Q) over the random magnetic field distribution [j9|]T7|] . The magnetic field in different 
domains is randomly oriented with respect to the neutrino propagation direction. So the neutrino conversion 
probability depends on the root mean square (rms) value of the random magnetic field. With the use of 
the delta correlation for uncorrelated magnetic field domains of size Lq, the average of the random magnetic 
field is H&|0), 

(B ]l (t)) = (B ± (t))=(B l] (t)B ± (t))=0, (8) 
(B q {t)B n (h)) = (BpSijLoSit - h), (9) 



and 



{B i± (t)B j± (ti)) - (Bl)S tJ L 5(t - h). (10) 



The rms value of the averaged magnetic field is given as B rrns = \/ (B 2 ). Let us assume, 

• H a b{t) is sufficiently close to iJ a b(^2) so that 

/ f 1 H ab {t) Hd{ ^ )H ^ 2) P^dhdh ~ /" r H d {t 1 )H d (t 2 )P{t)dt 2 dt 1 , (11) 
Jo Jo n a b[t2) Jq Jq 

and 

• P(t) has no correlation with H d {t) and H ab {t). 

Thus averaging both sides of Eq.(|j) and using the above magnetic field correlations we obtain 

V(t)=-2H 2 ab (0)([ (2P(t 1 )-l)dt 1 ) -T ± (2P(t)~l))~H 2 (0)(f P^dh) - 2T\\ ( f Pfodtt). (12) 
Jo Jq Jo 
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where (P) = V and the quantities T± and Ty are the longitudinal and transverse damping parameters given 

by 



4 9, „ , „ 1 

2 r>2' 



r_L = -^(B 2 )L , and r,| = -fil ff (B 2 )L . (13) 

Differentiating once again to Eq. (|l2|) with respect to t we obtain the following second order differential 
equation 



where the quantity 



V(t) + 2T T V(t) + u 2 s V{t) - 2H 2 ab (0) = 0, (14) 



u 2 s = 4H 2 ab (0) + H 2 (0), (15) 



is square of the spin rotation frequency of the neutrino in the medium, = T± + Tu and it satisfies the 
boundary conditions "P(O) = and "P(O) = T±. The solution to this differential equation is given by, 



The Eq. (|l^) is the general solution for the average conversion probability having random fluctuation in both 
transverse and longitudinal mode of the magnetic field |l5| ] . For no random fluctuation in the magnetic field 
(r_i_ = r ii = 0) we get back the standard MSW type solution for the neutrino conversion probability 

pu^JSA^A (17) 

Let us consider the propagation of a system of active (doublet) and light sterile (singlet) neutrinos (y e — > v s ), 
with masses mi and ni2, mixing angle 9, and no transition magnetic moments, in the presence of a magnetised 
plasma. Then the Hamiltonian in the evaluation equation Eq.(^J) will be 

'V- Acos26» + ^ e// k.B/fc Asin26>/2 S 
A sin 26»/2 



(18) 



where A = (m 2 — m\)/2E. For active-sterile neutrino conversion the resultant vector potential experienced 
by v e is given by 

V = V2G F n e {3Y e + AY Ve - 1), (19) 

where Gp is the Fermi coupling constant, n e is the electron density in the medium and Y e and Y Va are the 
electron and v e abundances respectively in the medium. For neutrino propagating along the z axis the V ax i a i 
is 

k z 

Vaxiai = HeffB z — - (20) 



The quantity Heff for v e - > Vs is given by 



eG F P F 

^ S = 7f^' (21) 
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and Pp is the Fermi momentum of electron. As we are considering the neutrino magnetic moment /transition 
magnetic moment to be zero, the perpendicular component of the damping term T± will vanish. Thus putting 
r± = in Eq.©, the average conversion probability for v a — > v s will be 



A 2 sin 2 28 / r ,, t L sinhUT 2 - uj t] 



V(t) = ( 1 - <| r„ j=L= + cosh [v /r 2 - J , } 1 . (22, 

where 

ujj = (Hj(0) + A 2 sin 2 26) (23) 

is square of the flavor conversion frequency of the neutrino in the medium. In a supernova medium with a 
strong random magnetic field satisfying the condition T|| 3> w/ the average conversion probability in Eq.(p2|) 
will be approximately 

A 2 sin 2 28 / ,2. , or \ , 

V{t) ~ 2uj 2 V ~ 6 ) ' (24) 

and for weak field limit i.e. T\\ <uj it will be 

A 2 sin 2 20 , r , s 

P(<) ~ 2 ^ 2 (1 ~ e_ C0SW/ ^ • (25) 

These are shown previously by Semikoz and Valle in ref |l3j . Using the strong field limit for the active-sterile 
neutrino conversion in a supernova, the limit on Am 2 and sin 2 28 are obtained from the supernova cooling 
llfj and by thermalisation of the sterile neutrinos in the early universe hot plasma JL3[| . 



B. method - II 

In the previous section we have derived the master equation for the average conversion probability of 
neutrino in the medium with a magnetic field. But in that calculation we have neglected the correlation 
of P(t) with Hd(t) and H a b(t). In the previous section we have only taken the average of the equation 
for P(t) and not for R(t) and I(t) to derive the master equation for the conversion probability. For the 
neutrino conversion length greater than the domain size (l con v >> Lq), a neutrino will cross many magnetic 
field domains before it flips its helicity. Thus the neutrino will experience an average field before it flips 
its helicity. So one can average the propagation equations (^|) and (^) over the random magnetic field 
distribution Let us define the average of the functions (P(t)) = V{i), (R(t)) = K(t) and (I(t)) = l(t). 

Because of the averaging the average probability V(t) will only depend on the even powers of the magnetic 
field correlation. Using the average functions in Eqs.(|) to @ we obtain 

±(t) = (H ab (t)(2P(t) - I)) + (H d (t)R(t)) 

= H ab (0)((2P(t) - I)) + (H ab (2P(t) - I)) 

+H aa {0)(R(t)) + (H aa (t)R(t)), (26) 

and 
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K(t) = -H aa (0)(I(t)) - (H aa (t)I(t)) (27) 
respectively. Using the delta correlations for the magnetic fields as shown in Eqs.(||) to ( |To| ) we obtain 

(H' d (t)R(t)) ~ -2r n X(t); (H' ab (t)P(t)) ~ -r x I(t), (28) 

and 

(ff^(f)J(t)>^(2P(t)-l). (29) 
Putting these values in Eqs.(p6|) and ( p7|) we obtain 

±(t) = H d (0)K(t) + H ab (0)(2T(t) - 1) - 2(r_L + r,|)Z(t). (30) 

and 

U(t) = -H d (0)X(t) - 2T ll K{t). (31) 
Using these equations in Eq.(^) we obtain 

T{t) = -2H ab (0)l(t) - T ± (2V(t) - 1) . (32) 

For convenience let us define 

J(t) =e- 2 ( r »+ ri)t X 1 (i), (33) 

and 

H(t) = e- 2r n*7£i(i), (34) 
Using Eqs.(||) and @ in Eqs.® to (||) we obtain 

±i(t) - ^(Oje^'Rjfi) + ff a6 (0) e 2 ( r -+ r ii> t (2P(i) - 1), (35) 

^ 1 (i) = -i7 d (0)e- 2r ^X 1 (t), (36) 

and 

P(t) = -2i7 Qb (0)e- 2 ( r ^ +r ii)%(i) - T x [2V{t) - 1) . (37) 



Differentiating Eq. (|37|) twice with respect to t and putting the value of X\ (t) from Eq. (|35|) we obtain the 
master equation for the average conversion probability as, 

v{t) + 4 (r_L + r,|) v(t) + 4 ^3r ± r,| + ri + rjj + ^\ v{t) 
+ 8 (rir„ + r x r| + a* (o)r,| + m 
- 4 ( rir,, + r x rf + # 2 6 (o)r„ + SjEi) = , (as) 



with the boundary conditions 7>(0) = 0, V(0) = Yj_ and P(0) = 2H% b (0) - 2Y\. By switching off the 
damping terms in Eq.(38), we will give get back the same solution for P(t) as shown in Eq. (|l7|) . The 
solution to the above third order differential equation is given in Eq.(25) of the ref and it is of the form 

= -+!/(*). (39) 

Henceforth for further reference we will refer the Eq.(25) of ref (T^|. The solution obtained is very complicated 
and it is difficult to conclude any thing from the solution. On the other hand the interesting part of the 
solution is that, the positive definiteness of the average conversion probability (0 < V < 1) requires the 
following condition to satisfy 

^ > ^(ri + r|-r x rn), (40) 

irrespective of the form of neutrino potential and the magnetic field, which essentially shows that neutrino 
crosses many domains leads to this requirement. Considering the magnetic field in the early universe hot 
plasma and the core of the newly born neutron stars to be purely random in nature, we studied the conversion 
v e L — > v e R in these medium in a previous paper |^9|. 



III. COMPARISON OF BOTH THE METHODS 

Let us consider the active sterile neutrino conversion in the magnetised medium, with neutrino mixing 
and fx = (r_i_ = 0) as shown in (p"8|). Only the fluctuation in the parallel component of the magnetic field 
will contribute. So for this case the condition in Eq.(|40|) will be modified to 

coj > -JL. (41) 

This implies u>f > Y\\. The average conversion probability in Eq.(p4|) for strong random magnetic field limit 
is derived with the approximation Y\\ 3> w/. On the other hand we obtain from Eq.(f4l|), completely the 
opposite one. Also the later condition further shows that, what ever may be the strength of the random 
component of the magnetic field Ty, it should always satisfy the condition in Eq. (|4l|) . Apart from that, 
this condition comes automatically from the positive definiteness of the probability. This opposite situation 
arises because of the averaging procedure. In the first case we only take the average of the Eq.(||) and use 
the functions I{t) and R(t) in that. On the other hand in the second method, we average the equations for 
P(t), I(t) and R(t) separately, which are shown explicitly. In the first method it is assumed that H a b(t) is 
sufficiently close to ffa&fe) so that in Eq. (|TT1) , they cancel. But in general this is not true for arbitrary t 
and t2- Such terms are not there in the second method to derive the average probability equation. 

Comparison of both the master equations Eq.(jlJ) and Eq.(^8|) shows that, in the first there is no mixing 
between the longitudinal rn and transverse Y j_ damping terms. Also the damping terms do not mix with 
the background terms. On the other hand, there is mixing between the damping and background terms, and 
the probability equation becomes more complicated. 
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IV. RESULTS AND DISCUSSIONS 



Let us estimate the range of Am 2 and sin 2 29 from the the above inequality in Eq.([4l|) in a supernova 
environment. Inside the supernova the vector potential for the process v e — > v s is 

V ~ 4 x 10 _ Vi4(3Y; + 4Y„ e - l)MeV ~ 4.48 eV, (42) 

where Ye ~ 0.3, Y Ve ~ 0.06 and the quantity pi4 ~ 8 is the density in the supernova core in units of 
10 14 gm/cm 3 . Now let us assume that the random fluctuation in the magnetic field is much large compared 
to the constant background part i.e. ^Bn <C T\\ as considered in ref [jl0|. The flavor precession frequency 
for u e — > v s is given by 



"t = \W-%£<*»W 2 +(^-) sin 2 20. (43) 



2E ' \ 2E 

Inside the supernova core, neutrinos have energy in the range 30 to 100 MeV. Then considering E ~ 100 
MeV, we obtain 



/ Am 2 / Am 2 \ 2 

uj f = W(4.48- 0.5 x 10- 8 — ^-cos26») 2 + 0.25 x 1Q- 16 ( sin 2 26 eV. (44) 



Inside the supernova core for u e 



Me// ~ 4.3 x 10- 13 MS ) • 14.-.) 



where fiB is the Bohr magneton. The quantity /i e // has the same dimension that of the magnetic moment, 
but it has nothing to do with the magnetic moment as it does not change the helicity of the particle. The 
Fermi momentum of the electron inside the core is 

P F ~ 320(y e pi 4 ) 1/3 MeV ~ 428 MeV. (46) 

Thomson and Dunkan have argued that very strong magnetic fields might be generated inside the supernova 
core due to small scale dynamo mechanism. If these fields are generated after core collapse, then it could be 
viewed as random superposition of many small dipoles of size Lq ~ 1 Km. Then the longitudinal damping 
parameter is simplified to 

r,| ~ 9.6£ 2 4 eV, (47) 

where B rms is expressed in units of 10 14 Gauss. Putting LOf and T\\ in Eq. (pl|) , we can find the ranges 
of Am 2 and sin 2 29 for which the condition in Eq.(f4l|) is satisfied. These are shown in the contour plots 
for different values of the magnetic fields in Figure 1. Figure 1 (a) shows that, a very narrow range of 
Am 2 /eV 2 (3.16 x 10 s — 1.6 x 10 9 ) is excluded for <J sm 2 26>0.63 x 10~ 3 for the random magnetic field 
strength B rms = 0.1 x 10 14 Gauss. By increasing the strength of the random magnetic field to 0.5 x 10 14 
and 0.63 x 10 14 Gauss, we see in Figure 1 (b) and (c) that the width of the excluded range of Am 2 /eV 2 has 
increased. The spread of Am 2 /eV 2 towards smaller values is large compared to the large values. Also the 
large values of sin 2 29 are excluded. In Figure 1 (b) we see that 3.5 x 10 s <, Am 2 /eV 2 <, 1.6 x 10 9 is excluded 
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and for this the excluded range of sin 2 29 is <S sin 2 28 <J 0.4 (for B rms = 0.5 x 10 14 Gauss). In Figure 1 (c) 
it is shown that for B rms — 0.63 x 10 14 Gauss, the excluded region of Am 2 /eV 2 is 1.6 x 10 7 — 1.6 x 10 9 and 
in this range all the values of sin 2 28 are excluded. It shows that, two distinct allowed regions for Am 2 /eV 2 
are there in both sides of the excluded region. In these allowed regions all ranges of sin 2 28 are allowed. But 
so far as Am 2 value is concerned, it is interesting to consider only the region left to the excluded curve, 
because it corresponds to smaller Am 2 /eV 2 values (Am 2 <, (keV) 2 ). Going from Figure 1 (c) to (d) (for 
B r ms = 10 14 Gauss) we observe that, the left arm of the curve vanishes and the right arm spreads towards 
the higher values of Am 2 /eV 2 for all ranges of sin 2 29. This implies that small values of Am 2 /eV 2 are 
excluded. By further increasing the strength of the random magnetic field, we have shown in Figure 1 (e) 
that, (for B rms = 10 16 Gauss) all values of Am 2 /eV 2 1.6 x 10 13 are excluded for all ranges of sin 2 29 
values. This range of Am 2 /eV 2 are obviously not at all interesting, because it corresponds to very high 
values of Am 2 /eV 2 , and no astrophysical, cosmological and laboratory observations favour this. Also if 
we consider small values of magnetic field (B rms < 10 13 Gauss) then we found that, the excluded region 
vanishes and all the parameter ranges are allowed. Thus if we consider only the effect of random magnetic 
field in a magnetised electron plasma along the neutrino propagation direction, then magnetic field should 
satisfy B rms <, 10 14 Gauss so that small Am 2 /eV 2 ranges should not be excluded. On the other hand for 
Brms < 10 13 Gauss we found that, all the parameter ranges are allowed. As we have shown in Figure 1 
(a), only a very narrow strip of Am 2 /eV 2 values are excluded for <, sin 2 29 < 0.4 and this narrow strip 
vanishes for smaller values of B rms . For Am 2 and sin 2 29 very small, we can see in Eq.(^) that ujf ~ 2.12 
eV and the condition Eq.([Il]) will give B„ns < 0.4 x 10 14 Gauss. Thus we found that for magnetic field 
Brms Ss -63 x 10 14 Gauss, all the interesting ranges of parameters are allowed. Also this implies that the 
maximum value of Am 2 can be in the (keV) 2 range. 



V. CONCLUSION 



For neutrino propagating in the magnetised plasma of the supernova, will experience an axial potential 
which is proportional to the scalar product of the neutrino momentum and the magnetic field. We derived the 
average conversion probability for active neutrino goes to active or sterile one by using a method developed by 
Scmikoz ct. al, and another by one of us where it is assumed that the magnetic field has a random fluctuation 
over the constant background. In the second method for neutrino having non-zero mixing and zero magnetic 
moment we found the condition for the positive definiteness of the average conversion probability for the 
process v e — > v s . Using this condition as the basis, we found the excluded/included ranges of Am 2 and 
sin 2 29 for different values of the magnetic field. Our calculation shows that for neutrino to have maximum 
mass in the keV range or less, the random magnetic field along the neutrino propagation direction should 
not be very large (B rm s 5s 0.63 x 10 14 Gauss). Because it is observed that for large magnetic field, small 
values of Am 2 are excluded. 



9 



(a) 



-io L 



8.5 8.6 8.7 8.8 8.9 9 



(b) 



-io L 



7 7.5 8 8.5 9 9.5 10 



(c) 



7 7.5 8 8.5 9 9.5 10 



(d) 



-10 



-10 



10 



(e) 

or ' "~ 



-2 



-4 
-6 



-8 



-io L , , , j 

10 11 12 13 14 

x 

Figure 1: We have denned x — log[Am 2 / eV 2 ] and y = log[sin 2 29].(a) It is for B rms — 10 13 Gauss and 
the shaded region is excluded, (b) for B rms = 0.5 x 10 14 Gauss and the region inside the curve is excluded, 
(c) B rms — 0.63 x 10 14 Gauss and the region inside the curve is excluded, (d) for B rms = 10 14 Gauss and 
left side of the curve is excluded and (e) B rms = 10 16 Gauss and the left side of the curve is excluded. 
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